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Abstract 

(N . 

It is shown that non-helical (more precisely, parity-invariant) flows 
. capable of sustaining a large-scale dynamo by the negative magnetic 

eddy diffusivity effect are quite common. This conclusion is based on 
numerical examination of a large number of randomly selected flows. 
£h ■ Few outliers with strongly negative eddy diffusivities are also found, 

j> | and they are interpreted in terms of the closeness of the control param- 

eter to a critical value for generation of a small-scale magnetic field. 
Furthermore, it is shown that, for parity-invariant flows, a moderate 
separation of scales between the basic flow and the magnetic field of- 
ten significantly reduces the critical magnetic Reynolds number for 
the onset of dynamo action. 
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1 Introduction 



The kinematic dynamo effect is generally believed to be facilitated by a 
separation of spatial scales between the flow and the magnetic field. This 
idea has its roots in the discovery of the a-effect (Steenbeck et ai, 1966; 
Krause and Radler, 1980) for non-parity- invariant flows (an instance of which 
are helical flows). Indeed, recent results indicate that helicity and scale 
separation can lead to experimental dynamos (Gailitis et ai, 2000; Busse, 
2000, Stieglitz and Muller, 2001). 

Lanotte et ai (1999) showed that scale separation can also be useful for 
flows possessing parity invariance and thus lacking a-effect. They found 
a first instance of a three-dimensional parity-invariant flow, the "modified 
Taylor-Green flow", which can amplify a large-scale magnetic field by a 
negative eddy diffusivity mechanism, as conjectured by Kraichnan (1976). 
Examples of negative eddy diffusivities when the magnetic field depends 
on all three coordinates, while the flow depends only on two (two and a 
half dimensional dynamos) were actually known since the work of Roberts 
(1972). New instances of large-scale dynamos with the Roberts flow were 
found more recently by Tilgner (1997) and Plunian and Radler (2000, 2001). 
Emergence of the ct-effect and the negative eddy diffusivity effect can be 
demonstrated by asymptotic theory involving, in principle, very large scale 
separation. Both the aforementioned experiments and numerical simulations 
with helical (Galanti et al., 1992) flows indicate that even a moderate scale 
separation of a factor two can help for obtaining a dynamo (a physical in- 
terpretation for the case of the ABC flows has been proposed by Archontis 
(2000)). 

The following argument due to Fauve (1999) shows, how crucial for ex- 
perimental dynamos it is to lower the critical magnetic Reynolds number 
R m = LV/t]. Here, L is the overall scale of the apparatus, V is the typical 
velocity and rj is the (molecular) magnetic diffusivity. Given the typically 
very small value of the magnetic Prandtl number, the flow is expected to be 
strongly turbulent. Hence, the kinetic energy dissipated per unit time in the 
apparatus is P ~ p(V 3 /£)L 3 , where p is the fluid density and I is the typical 
spatial scale of the velocity. Thus, P ~ prfR^KeL), where e = £/L is the 
scale-separation factor. Therefore, for a given fluid (rj), geometry (e) and 
size (L), the required pump power behaves as the cube of the attainable R m . 
It is generally believed that the cost of an experiment scales roughly as its 
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volume, i.e. ~ L 3 ; thus, if there are practical limitations on the maximum 
pump power available, the cost will scale as the ninth power of the attainable 
R m . All this makes it desirable to design apparatus with the lowest possible 
critical magnetic Reynolds number. 

Our intention is to begin a systematic investigation of the kinematic dy- 
namo effect for parity-invariant flows with large or moderate scale separation 
with emphasis on (i) the genericity of the phenomenon of negative eddy diffu- 
sivity and (ii) the lowering of the critical magnetic Reynolds number. Note, 
that parity-invariant flows are at least as "natural" as helical ones. Indeed, 
the Navier-Stokes equations are parity- invariant. Hence, even if this par- 
ticular symmetry, like most other symmetries, is broken by the mechanisms 
generating the flow, we can expect the symmetry to be dynamically restored 
in a statistical sense at high Reynolds numbers (Frisch, 1995). Note also that 
the existence of three-dimesnional flow with negative eddy diffusivity is a re- 
cent discovery; contrary to the alpha effect which can be derived in purely 
analytic fashion by working, e.g., at low magnetic Reynolds numbers or with 
a delta-correlated (in time) velocity, only (very CPU-demanding) numerical 
solutions of the equations arising from multi-scale techniques can demon- 
strate the existence of negative eddy-diffusivities and can show that it is not 
a pathologically unusual effect. We shall work here with time- independent 
flows, which we may think of as being averages of turbulent flows. We are of 
course aware that time dependence of the flow and, generally more realistic 
choices of the basic flow could somewhat modify our conclusions; such issues 
are left for future work. 

The paper is organized as follows. 

In Sections |^ and |3| we consider the case where the spatial scale of mag- 
netic field is much larger than that of the flow. Then e is a natural small 
parameter, enabling one to apply asymptotic methods. Multiscale expansions 
in hydrodynamics may be viewed as adaptations of general homogenization 
methods used to determine effective large-scale transport coefficients of pe- 
riodic or random materials (Bensoussan et al, 1978; Kozlov, 1978; Zhikov et 
al, 1979). Usually the effective coefficients cannot be obtained without solv- 
ing suitable auxiliary small-scale problems. In three space dimensions this 
may require significant computational resources. In Lanotte et al. (1999) 
equations for the large-scale component of magnetic fields are derived along 
the same lines as in previous works on negative eddy viscosity in ordinary in- 
compressible hydrodynamics (Dubrulle and Frisch, 1991; Gama et al, 1994; 
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Wirth et ai, 1995). Here we use a similar but a distinct approach based on 
the work of Vishik (1986, 1987), which yields a complete decomposition of 
magnetic modes and their eigenvalues (in Appendix [A] arbitrarily high-order 
terms of such expansions are derived). 

To investigate, how commonly the phenomenon of negative magnetic eddy 
diffusivity is exhibited by fully three-dimensional flows, we follow in Section |3] 
a strategy similar to the one, which Gama et al. (1994) applied to estimate 
the abundance of two-dimensional flows with negative eddy viscosity. By 
solving numerically the relevant auxiliary problems we determine the mini- 
mum over directions of the (anisotropic) magnetic eddy diffusivities for two 
ensembles of flows with random Fourier components. This way we determine 
the fraction of flows such that, when decreasing molecular magnetic diffusiv- 
ity, a large-scale negative eddy diffusivity instability appears before the flow 
becomes a small-scale dynamo. An interpretation of "outliers" with strongly 
negative eddy diffusivity, found in computations, is given in Appendix |B[ 
Since it is important for our purpose to analyze as many instances as possi- 
ble, we have used a new efficient algorithm for solving the elliptic auxiliary 
problems arising from the two-scale expansion. This method is an outgrowth 
of the Chebyshev iterative method developed by Podvigina and Zheligovsky 
(1997) for self-adjoint negative definite linear operators. The operators aris- 
ing in hydrodynamical problems are not self-adjoint, but they naturally split 
into a dominant self-adjoint negative definite part (the molecular diffusion 
term) and the remaining transport term. We can then combine Chebyshev 
methods with operator splitting to reduce the harm stemming from com- 
plex eigenvalues. This method and its theoretical foundations are presented 
in detail in Zheligovsky (2001), available only in Russian. For the reader's 
convenience, the essentials are given in Appendix [Cj 

In Section |] we turn to cases with moderate scale separation (i.e. when e 
is finite). First, we consider instances of flows for which we already know that 
there is a negative magnetic eddy diffusivity, and we determine windows in e 
over which dynamo action is present. Second, we examine the case e = 1/2, 
that is, we assume that the magnetic field has a spatial period twice that 
of the flow, and study how the critical magnetic Reynolds number R m is 
affected, compared to the case with no scale separation (e = 1). In both 
types of computations we exploit the fact that, for space-periodic boundary 
conditions, eigenmodes have the form 

H = e^ x h(x;e,q), (1) 
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where q is a prescribed wavevector. This enables one to solve a problem 
involving only small scales, but parameterized by e and q (cf. Roberts, 1972; 
Frisch et ai, 1986). Solution of the eigenvalue problem is more CPU- intensive 
than the solution of auxiliary problems involved in the asymptotic case (be- 
cause efficient preconditioning techniques are available for the latter). Nu- 
merical methods and codes used for the eigenvalue computations are adapted 
from Zheligovsky (1993). 

In Section [| we summarize our main results. 

2 Asymptotic expansion of magnetic modes 
and of their eigenvalues 

Here we briefly review the asymptotic expansion arising in the kinematic 
dynamo problem, when the scale ratio e is small. The conducting fluid flow 
v(x) is assumed to be time- independent, and the problem therefore reduces 
to the eigenvalue problem for the magnetic induction operator 

r?V 2 H + V x (v x H) = AH (2) 

in the space of solenoidal vector fields: 

V ■ H = 0. (3) 

Here rj > is magnetic molecular diffusivity, H(x) is a magnetic mode and 
A is the associated eigenvalue (the real part of which is here referred to as 
the growth rate of the mode; negative growth rates correspond to decaying 
modes). The flow is assumed incompressible: 

V ■ v = 0, (4) 
27r-periodic in space and parity-invariant: 

v(x) = -v(-x). (5) 

The magnetic mode H(x, y) is supposed to depend on the fast variable 
x and on the slow variable y = ex, and is 27r-periodic in both fast and slow 
variables (the velocity v is independent of slow variables). Introduction of the 
two spatial scales results in a "decomposition of the gradient" V — > V x + eV y 



5 



in (0) and ([3]). Here and in what follows the subscripts x and y refer to 
differential operators in fast and slow variables, respectively. 

Let (•) and {•} denote the mean (over the periodicity cube) and the 
fluctuating part of a vector field, respectively: 

(f) = (2vr)- 3 / f(x,y)dx, {f} = f-<f>. (6) 

Parity invariance (|5|) implies (v) = 0. 

Solutions to the eigenproblem (H)-© are sought in the form of asymptotic 
series 

oo 

H = $:h ri (x,y)6 n , (7) 

n=0 

oo 

A = £ \ n e n . (8) 

n=0 

As shown in Appendix [A], substitution of these expansions in (|^) leads to a 
hierarchy of equations. By application of standard techniques for singular 
perturbation expansions, terms in the expansion eigenvalues are successively 
determined from solvability conditions, obtained here by just averaging the 
equations. In particular, it can be shown that Ao = Ai = 0, and the leading 
term A2 in the series (|8]) can be determined from the following eigenvalue 
problem 

3 3 ft/h k \ 
77Vj<h > + V y x £ £ B mk ^pl = A 2 (h ), (9) 

k=l m=l °Vm 

subject to the solenoidality condition 

Vy " (h ) = 0. (10) 

Thus, the leading terms of the expansions in powers of e of the mean mag- 
netic field and of the eigenvalue turn out to be, respectively, an eigenfunc- 
tion and the associated eigenvalue of a partial differential operator comprised 
only of second order derivatives, which may be regarded as representing an 
anisotropic diffusion. However, this operator is not necessarily dissipative: it 
may have eigenvalues with positive real parts. In this context one then may 
speak of negative eddy diffusivity. 

As shown in Appendix |A| eqn. ( |A.13 ), the coefficients D m fc of the ho- 
mogenised equation (|9|) are given by 

D mfc = (v x r mk ), (11) 



6 



where the tensor field T mk is obtained by solving a set of two auxiliary prob- 
lems: 

,V 2 S t + Vx(vxS t ) = -- (12) 



dx k 



and 



AC 

V V 2 T mk + V x (v x T mk ) = -277—^ + v m (S k + e k ) - vS™ (13) 

ox m 

which will be referred to as the first and the second auxiliary problem, respec- 
tively. Here, is the unit vector in the fc-th cartesian coordinate direction 
and superscripts denote cartesian components of vectors. For solutions to 
( |12"D and ([H|) to exist, we require that zero should not be an eigenvalue of 
the magnetic induction operator 

£H = i)V 2 H|Vx (v x H), 

restricted to small-scale solenoidal vector fields with a vanishing mean; by 
"small-scale" we mean independent of the slow variable. Generically, this 
condition is satisfied. 

Problems (|T2|) and ([13|) are formulated with a 27r-periodicity bound- 
ary condition for the vector fields Sfe(x) and T mk (x.), which should also 
have a vanishing mean. It follows that the S k s are parity-antiinvariant 
(Sfe(x) = Sfc(— x) ) and solenoidal (V-S& = 0), and that the T mk s are parity- 
invariant (r mfc (x) = -r mfc (-x) ). 



Since the homogenized operator defined by fl5|)-(|TUD has constant coeffi- 
cients, its eigenvectors are Fourier harmonics: (ho) = he* qy , where h and 
q are constant vectors. The wavevector q is taken real because we are only 
interested in modes uniformly bounded in the entire space. Substitution in 

(D-n y ields 

3 3 

7?|q| 2 h + qx^ V mk h k q m = -A 2 h, (14) 

fc=l m=l 

h-q = 0. (15) 
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3 Abundance of flows producing negative mag- 
netic eddy diffusivity 

Following Gama et al. (1994), we now ask how common is the phenomenon 
of negative diffusivity. 

Time-independent periodic flows, satisfying ([§)-(|5]), are represented here 
by their Fourier series 

v = E v ke ikx (16) 
|k|=i 

with random amplitude harmonics. For the velocity to be real, Vk = v_k, 
and parity invariance (|5|) requires Revk=0. 

We study two ensembles of flows with different types of fall-offs of the 
energy spectrum Ek- The latter is defined as the total energy contained 
in harmonics, whose wavevectors k belong to the K-th spherical shell, i.e. 
K — 1 < |k| < K, where K > is integer. When Ek oc e~^ K , the flow is said 
to have exponential spectrum; when Ek oc K -1 , it is said to have hyperbolic 
spectrum. These have respectively rapid and slow fall-offs. The spectra are 
band-limited and set to zero beyond wavenumber N = 10 for the exponential 
case and beyond N = 7 for the hyperbolic case. The value £ = 1CT 2 / 3 is used 
throughout in the exponential case. 

The flows are produced by the following procedure (applied only to one 
half of the set of wavevectors, the second one being related by complex con- 
jugation). For k in the ball 1 < |k| < N, random imaginary vectors rk are 
generated with the imaginary part of each cartesian component uniformly 
distributed in the interval [—0.5, 0.5] (the RAN2 random number generator 
by Press et al, 1992, is used). Then the vectors rk are projected onto the 
plane perpendicular to k to gain incompressibility. Finally, the obtained vec- 
tors are normalized, so that the flow has the desired energy spectrum and the 
root-mean-square velocity is one, when the normalized vectors are employed 



as the Fourier coefficients Vk in (16). 

Each ensemble consists of 100 flows. For each flow, computations are 
carried out for three values of molecular diffusivity, rj = 0.3 , 0.2 and 0.1 . 
Note that when rj is too large, no dynamo action can take place at large or 
small scales because the eddy diffusivity will basically be equal to the molec- 
ular diffusitivity. (In contrast, an alpha-effect dynamo is still possible if the 
system is large enough). If rj is too small, small-scale dynamo instabilities 
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Table 1. Statistics of the eddy diffusivities for the considered ensembles of 
flows. 





Exponential spect 


rum flows 


Hyperbolic spectrum flows 




rj = 0.3 


r] = 0.2 


T] = 0.1 


rj = 0.3 


77 = 0.2 


r] = 0.1 


Re ^eddy < 


0% 


18% 


86% 


0% 


4% 


53% 


Re??eddy < T) 


83% 


96% 


98% 


30% 


63% 


94% 


complex i] eddy 


7% 


2% 


2% 


8% 


5% 


3% 



will have much larger growth rates than large-scale instabilities. In both 
extreme cases it is meaningless to study negative eddy diffusivity large-scale 
instabilities. The range of values of rj chosen, between 0.1 and 0.3, is an em- 
pirically determined tradeoff between these extremes. All computations are 
performed with the resolution of 64 3 Fourier harmonics, employing standard 
spectral methods based on fast Fourier transforms and dealiasing. (Each case 
requires from 10 to 30 minutes of CPU on one DEC Alpha EV6 processor, 
depending on the value of rj.) The power spectra of the solutions to fll~2|) and 
( |13|) turn out to decrease monotonically very quickly for K > 2. For r] — 0.1 
in the hyperbolic case (resp. exponential case), the decrease is typically by 
7-8 orders of magnitude (14-15 orders of magnitude) from maximum to min- 
imum. Note that our code reproduces the negative eddy diffusivity values 
obtained by an alternative numerical method by Lanotte et al. (1999) for the 
modified Taylor-Green flow. 

After the homogenized coefficients (|Tl~|) are determined, the minimum 
magnetic eddy diffusivity 

?7eddy = min(-A 2 (q) ) 
|q|=i 

is calculated as the minimum eigenvalue of the problem (p!4|) -(|T5|). (Note, 
that for a unit wavevector the eigenvalue is just the opposite of the eddy 
diffusivity.) For a given q this eigenproblem yields two eigenvalues satisfy- 
ing a quadratic equation; the eigenvalues can thus be complex conjugate (a 
similar phenomenon was observed for three-dimensional eddy viscosities; see 
Wirth et al, 1995). Then the minimum of the real part of the complex eddy 
diffusivity, which determines the maximum rate of instability of a large-scale 
magnetic perturbation, is sought. However, in our computations complex 
eddy diffusivities appear quite rarely. 
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(a) (b) 

Figure 1: Histograms of the minimum eddy diffusivity values r] eddy (horizontal 
axis) for ensembles of 100 parity-invariant flows with exponential (a) and 
hyperbolic (b) spectra. 
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Histograms of the obtained magnetic eddy diffusivities are shown in Fig- 
ure 0. The statistics of eddy diffusivities shown in Table 1 indicate that the 
lower the molecular diffusivity, the higher the chance that a flow will have a 
negative "eddy correction" (correction to the molecular value coming from 
the flow) and also the higher the chance that the resulting eddy viscosity will 
be negative. Note however that the eddy correction (given, at least implicitly, 
by (|TlD -(|l3D) has a complicated, generally non-monotonic, nonlinear depen- 
dence on the molecular diffusivity. Table 1 also shows that flows with the 
exponential energy spectrum are more efficient in lowering the minimum of 
the eddy diffusivity than those with the hyperbolic spectrum. This suggests 
that low wavenumber harmonics of the underlying flow are more important 
for lowering than high wavenumber harmonics. 

Examples of magnetic eigenmodes for 77 = 0.1 are shown in Figures |2] 
and |3|. They correspond to two samples having exponential velocity spectra. 
The upper of each figure shows regions where |v| 2 is less than 0.2% of its 
maximum, hence revealing the immediate neighborhood of stagnation points. 
Note that in consequence of parity invariance about the origin and of 2n- 
periodicity, within a basic 27r-periodicity cube there are in total eight points 
about which the flow is parity- invariant. These stagnation points are located 
at the vertices of a cube with side n and they form a uniform 7r-periodic grid in 
space. The magnetic structures seen in the lower left and right of each figure 
are clearly located near some of those stagnation points. These structures 
are magnetic flux ropes of the kind considered by Galloway and Zheligovsky 
(1994). Note that only the fluctuating part of the magnetic eigenmodes are 
shown (both the periodicity-box averaged magnetic field (ho) and the overall 
slow spatial factor e ieqx are omitted). In Figure |3] (ho) is very small (about 
two orders of magnitude smaller than the fluctuating field), whereas for the 
case of Figure ^| it is of the same order of magnitude as the fluctuating field. 

A physically important question is: which instability arises first, when 
magnetic molecular diffusivity is decreased, a large-scale or a small-scale one? 
Similarly to what was done by Gama et al. (1994), for each considered pair 
of a flow and a molecular diffusivity, it is checked that small-scale magnetic 
fields with a vanishing mean are stable, i.e. have negative growth rate (see 
Figures f| and ||) . 

Several "outliers" with relatively large negative eddy diffusivities can be 
seen on the histograms shown on Figure [L] for 77 = 0.1 . This can be ex- 
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Figure 2: First sample of velocity field and corresponding magnetic eigen- 
function shown as isosurfaces. Upper: |v| 2 at 0.2% of its maximum; lower 
left: |H| 2 at 60% of its maximum; lower right: |H| 2 at 40% of its maximum. 
Vertices of the bold solid line cube are stagnation points; the relevant ones 
are labelled a, b and c, the associated magnetic stuctures being labelled with 
the corresponding capital letters; primes and double-primes distinguish pairs 
of structures; subscripts such as 1 or 2 indicate the continuation of a given 
structure by periodicty. 19 



Figure 3: Second sample of velocity field and corresponding magnetic eigen- 
function shown as isosurfaces. Upper: |v| 2 at 0.2% of its maximum; lower 
left: |H| 2 at 20% of its maximum; lower right: |H| 2 at 8% of its maximum. 
Otherwise as in Figure 



13 



-0.4 -0.3 -0.2 -0.1 

Figure 4: Minimum eddy diffusivity redely (vertical axis) versus the growth 
rate of the corresponding dominant small-scale magnetic modes (horizontal 
axis) for rj = 0.3 (dots), 0.2 (pluses) and 0.1 (triangles) and 100 different 
flows with the exponential spectrum. 
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Figure 5: Same as Figure £| but with the hyperbolic spectrum. 
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Figure 6: Maximum growth rates (\ of small-scale magnetic modes (a), min- 
imum magnetic eddy diffusivity r/ e dd y (b) and product Ci^eddy (c) all plotted 
versus molecular diffusivity rj. Line width codes two flow realizations, corre- 
sponding to the two right-lower corner outliers in Figure U 
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plained as follows. We see on Figures f| and [5| that in the case of outliers 
the growth rate of small-scale magnetic modes (plotted horizontally) is nega- 
tive, but close to zero. Furthermore, we have checked that the corresponding 
eigenvalue is real. In other words, we are close to a non-oscillatory small- 
scale instability threshold. Now, we observe that the expression of the eddy 
diffusivity, derived in Section |2| and Appendix [A], involves the inverse of the 
small-scale magnetic induction operator. It is thus not surprising that large 
eddy diffusivities can be obtained. Depending on whether a certain inequality 
is satisfied at the threshold, this eddy diffusivity can be either large negative 
or large positive (see details in Appendix |FJ). We have checked numerically 
that the outlier eddy diffusivities are approximately inversely proportional 
to the small dominant eigenvalue (see Figure 

Note, that if we are exactly at the threshold, our two-scale expansion 
breaks down. A modified expansion can be carried out which takes into ac- 
count the enlargement of the kernel of the magnetic induction operator L. In 
the presence of a large-scale perturbation this can lead to an eigenvalue 0(e) 
rather than the 0(e 2 )-dependence obtained in the general parity-invariant 
case. In other words, the large-scale dynamics can be governed by a first- 
order differential operator. 

4 Generation of magnetic field with moder- 
ate scale separation 

For intermediate values of e, corresponding to moderate scale separation, 
direct evaluation of the most unstable magnetic modes is desirable. Compu- 
tations are simplified by the fact that eigenmodes have the form ([D), all the 
dependence on the slow variable y = ex being absorbed into an exponential 
factor e* eq ' x . Since the flow is assumed to depend only on the fast variable x, 
the eigenvalue problem is reduced to a small-scale one for a modified operator 
by changing V — > V x + ieq in the magnetic induction operator (^) and in 
the solenoidality condition @. 

To find the most unstable magnetic mode one must, in principle, allow for 
arbitrary orientations of the vector q. Throughout this section we assume 
only "binary" q's such that each component takes only the values or 1, 
so that e _1 may be regarded as the maximum scale-separation factor over 
cartesian directions. 
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-0.05 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 

Figure 7: Maximum growth rates of magnetic modes (vertical axis) versus 
the scale ratio e (horizontal axis) for different molecular diffusivities r] and 
one of four samples of a flow with the hyperbolic spectrum. 

Since our simulations are computationally demanding, we choose for each 
flow only one q-direction as follows. For a given flow and for the molecular 
diffusivity rj — 0.1, we obtain the most unstable q-directions in the asymptotic 
regime e — > from the results of Section It turns out that for the cases 
reported in this section, these directions are approximately parallel to binary 
vectors with one or two non-zero components. These binary vectors are then 
employed when e is finite for all values of rj studied. 

Figures show the dominant growth rates for four different flows se- 
lected out of the ensemble of flows with hyperbolic spectrum, which sustain 
negative eddy diffusivity. Actually, when q has integer components, it is only 
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0.1 0.2 0.3 0.4 0.5 

Figure 8: Same as Figure [7], another flow sample. 

necessary to make computations for < e < 1/2. Indeed, first note, that the 
modified magnetic induction operators for two values of e differing by 1 have 
the same spectrum and their eigenfunctions differ by a factor e* qx . Second, by 
parity invariance of the velocity fl5|), for two opposite values of e the spectra 
are also the same and eigenfunctions are related by parity H(x) — ► H(— x). 
(From this one easily deduces that the spectrum of the modified operator is 
complex conjugate, eigenvalues A and A being associated to eigenfunctions 
H(x) and H(— x), respectively.) Hence, outside the interval < e < 1/2 the 
plots can be continued by symmetry about the vertical line e = 1/2 and by 
1-periodicity. 

The results reported here are obtained with the resolution of 64 3 Fourier 
harmonics. For the "worst" case, r\ = 0.07, they differ only in the third 



19 



0.002 



0.001 



-0.001 - 



-0.002 - 



-0.003 



-0.004 




?7=0.069 



-" 77=0.070 



=0.071 



=0.072 



=0.073 



7?=0.074 



77=0.075 



0.1 0.2 0.3 0.4 

Figure 9: Same as Figure |7|, yet another flow sample. 



significant digit from those obtained with the resolution 32 3 . Note, that 
energy spectra of the computed magnetic modes decrease by at least 4 orders 
of magnitude for r\ = 0.13 and by at least 5 for r\ = 0.3 , indicating that the 
employed resolution is adequate. (For rj = 0.07 the decrease is only by 3 
orders of magnitude, which may seem marginal. However, it is well-known 
that eigenvalues are less sensitive to insufficient resolution, than the fine 
structure of eigenmodes.) It is clear from Figures |^-|10| that for all instances 
where the eddy diffusivity is negative (i.e. the graph go to zero by negative 
values as e — > 0) there is also an instability for any finite e. The converse 



is, however, not always true; for example, for the last flow (Figure |10|) and 
T] = 0.10 there is a dynamo for 0.26 < e < 0.5, which disappears below 
T] = 0.26 (the eddy diffusivity is positive). 

Furthermore, inspection of Figures |3|-|T0| reveals that there can be one or 
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Figure 10: Same as Figure |7|, the last flow sample. 

more windows in scale separation, where the dynamo (positive growth rate) 
is present, whereas in Figures [7] the dynamo appears to be present for all 



scale separations, provided the eddy diffusivity is negative. In Figure [10] the 
two windows belong to two different analytical branches of eigenvalues and 
eigenmodes, with the transition occurring near e = 0.3 . 

Note, that all the growth rates vanish (quadratically) when e — > 0. The 
reason is that the small-e eigenvalues, obtained by the two-scale expansion 
of Section 0, are analytic continuations of the zero eigenvalue in the small- 
scale problem. The presence of the eigenvalue zero is due to the structure of 
the magnetic induction operator, which is the curl of a first-order operator. 
(Hence, its adjoint, acting also in the space of divergenceless functions 

C* : Hh r/V x H + (V x x H) x v - V x p 
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Figure 11: Critical magnetic Reynolds numbers for the onset of instabil- 
ity of small-scale magnetic modes (horizontal axis) versus critical magnetic 
Reynolds numbers for magnetic modes of twice the spatial period (vertical 
axis). 28 flows with the hyperbolic spectrum are shown. The dashed line 
corresponds to equal Reynolds numbers. 



has constant vectors in its kernel.) Conversely, it can be demonstrated by 
spatial averaging of £H = AH, that for a mode with a non-vanishing mean 
the associated eigenvalue equals 0. (Note, that in Section we restricted 
the operator to functions of vanishing spatial mean and thus the growth rate 
zero was absent from Figures |4] and [|.) 

We now turn to the case, where there is a finite scale separation of a 
factor two (e = 1/2). The question we address is: how does a mere doubling 
of the spatial period of magnetic modes (relative to that of the basic flow) 
affect generation? More precisely: to lower the critical magnetic Reynolds 
number, does it help to introduce such a scale separation? We systematically 
compare critical numbers when e = 1 (no scale separation) and e = 1/2 for an 
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ensemble of flows with hyperbolic energy spectrum. We define our magnetic 
Reynolds numbers R m = LV/rj in terms of the largest scale available L = e~ l 
(the inverse of the minimum wavenumber of the magnetic field) and of the 
r.m.s. velocity V. In computations V = 1 and hence R m = (er?) -1 . 

When e = 1/2 we face again the problem, mentioned earlier, of choosing 
the vectors q. For each flow we now search for the critical R m , allowing all 
7 possible binary q ^ 0. More precisely, we find an approximate (minimum) 
critical R m and the corresponding q by running our code with 32 3 Fourier 
harmonics, and then refine the critical R m with 64 3 Fourier harmonics for 
the obtained q. 

The secant method is used to find stationary modes. Iterations arere 
terminated when the real part of the dominant eigenvalue was below 10~ 5 in 
magnitude. This usually requires 4-5 iterations. 

As far as the accuracy of the calculation is concerned, note that in all 
e = 1/2 runs with R m > 25 the energy spectrum decreases by at least 
4 orders of magnitude, an indication that 64 3 harmonics provide sufficient 
resolution. In the most demanding runs, with e = 1 and R m > 50, the 
energy spectra decrease only by 2-3 orders of magnitude but, as pointed 
out before, eigenvalues are computed nevertheless quite accurately (critical 
i? m 's, computed with the 32 3 and 64 3 resolutions, differ in the third significant 
digit). 



Figure |TT| shows all the critical -R m 's calculated for 28 cases. Two cases 
not shown have critical R m = 30.7 and 25.3 when e = 1/2, but do not display 
dynamo action for any R m < 100 when e = 1. The results demonstrate that 
scale separation lowers the critical R m in all but 9 cases out of the 30 studied. 
Note that Nore et al. (1997) have studied a flow (a variant of Taylor-Green) 
for which a scale separation e = 1/2 seems to lower significantly the critical 
magnetic Reynolds number, as is the case for the helical ABC flows studied 
by Galanti et al. (1992), Archontis (2000) and Dorch (2000). However the 
Nore et al. (1997) definition of the magnetic Reynolds number uses the scale 
of the flow and not that of the magnetic field. If instead our definition is 
used, their flow appears to be just a borderline case. 

5 Conclusion 

Three main results are presented in this work. First, in Section |3], we demon- 
strated that parity-invariant flows possessing a negative magnetic eddy diffu- 
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sivity are quite common. Specifically, we generated two random ensembles of 
flows; both had time-independent randomly chosen Fourier components, one 
with an exponentially decreasing energy spectrum, the other - with a power- 
law decrease (referred to as "hyperbolic" because the exponent is -1). One 
hundred independent selections were made in each case and three different 
molecular magnetic diffusivities were used, all corresponding to small-scale 
magnetic Reynolds numbers not exceeding ten. None of the 600 cases inves- 
tigated was a small-scale dynamo. For the lowest diffusivity considered, 86% 
of the "exponential" spectrum flows and 53% of the "hyperbolic" spectrum 
flows produced a negative eddy diffusivity and, hence, can act as a large-scale 
dynamo (for more details, see Table 1). 

Second, several instances with particularly negative eddy diffusivities were 
encountered. The theory of these "outliers" (given in details in Appendix ||) 
relates them to a near-critical state, where an eigenvalue of the small-scale 
dynamo problem is about to pass through zero (there is also an additional 
inequality ensuring that a large negative value is obtained). 

Finally, we presented numerical evidence that, for a parity-invariant flow, 
even a moderate scale separation between the flow and the magnetic field 
is generally beneficial for magnetic field generation. Specifically, we showed 
that the critical magnetic Reynolds number based on the largest scale present 
(that of the magnetic field) is often lowered when the scale of twice the spatial 
period of the velocity is introduced (70% of the instances studied had this 
property; see Figure |lTJ) . For a few flows we also studied the dependence 
on the scale-separation factor, when it is varied from a 1:2 ratio to a 1:20 
ratio. One instance (shown in Figure [10]) has a window lacking dynamo 
action between low and large values, where it is present. 

For helical flow, such as ABC, the lowering of the critical magnetic 
Reynolds number, when a moderate separation of scale is assumed, is un- 
derstood in terms of symmetry breaking, stagnation points and other ge- 
ometrical features (Archontis 2000). It would be of interest to find which 
geometrical features are important for the lowering mechanism when the 
flow is parity- invariant. 
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Appendices 



A Formal asymptotic expansion of magnetic 
modes and the associated eigenvalues 

We construct here complete asymptotic expansions of magnetic eigenmodes 
and of their eigenvalues for small values of the scale ratio e. Assumptions 
concerning the flow are listed in Section ^|. 

The kinematic dynamo problem reduces to the eigenvalue problem 

MH = AH. (A.l) 

Here H(x) is a magnetic mode, A is its eigenvalue, and 

MH = t?V 2 H + V x (v x H), 

acting in the space of solenoidal vector fields, is the magnetic induction op- 
erator (it acts like C, but in a larger domain). The magnetic mode H(x, y) is 
supposed to be solenoidal (^) and 27r-periodic in both fast and slow variables. 
Solutions to the eigenproblem ( |A.1| ) are sought in the form of asymptotic se- 
ries H = h n e n and A = A n e n - Denote 

H n = (h n ); G n = {h n } , 

where (■) and {•} refer to the mean and fluctuating part (see @ ), respec- 
tively. 

The expansion is constructed under the assumption that zero is not an 
eigenvalue of the operator C acting in the space of small-scale periodic 
solenoidal vector fields with vanishing mean. (By "small-scale" we under- 
stand "depending only on the fast variables".) 

The decomposition of the gradient V = V x +eV y , applied to the solenoidal- 
ity condition, expanded in powers of e and separated into mean and fluctu- 
ating part, gives 

V y • H„ = 0, (A.2) 
V x • G n + V y • G n _! = (A.3) 
for all n > 0. (H n = G n = is assumed for all n < 0.) 
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Substituting the expansions and the decomposition of the gradient into 
(|A.1|) , we obtain: 



E 

n=0 



CG n + V (2(V X ■ V y )G n _! + V y (H n _ 2 + G n _ 2 )) + V x x (v x H n ) 

n = 0. (A.4) 



+ V y x (v x (H 

n-1 + G„_i 

n—m 

(H m + G 

rn 

m=0 

We proceed by successively equating to zero the mean and the fluctuating 
part of each term of this series. In our context, averaging is equivalent to 
obtaining the solvability condition for a partial differential equation in fast 
variables. 

i. Transform the leading (n = 0) term of ( |A.4| ) using vector algebra 
identities and (W): 



£G + (H • V x )v = Ao(H + Go). 



(A.5) 



Averaging of this equation yields Ao = 0. Then the solution of (|A.5| ) can be 

expressed as 

G = £S fe (x)iJ fc (y), (A.6) 

k=l 

where the vector fields are the solutions of the first auxiliary problem (JT2|); 
the solutions exist since it is assumed that the kernel of £ is empty. (The 
first and second auxiliary equations were stated in Section 0.) 

By taking the divergence of the first auxiliary equation (|T2"D, one obtains 
V x • Sfc = 0. Since the velocity is parity-invariant, the domain of £ is a 
direct sum of two invariant subspaces, one of which is comprised of parity- 
invariant vector fields (i.e. f(x) = — f(— x) ), and the second - of the parity- 
antiinvariant ones (f (x) = f (— x) ). Thus, S is a parity-antiinvariant matrix: 

S fc (x) = S fc (-x). (A.7) 

ii. The e 1 term of ( |A.4j ), after the use of (|A.2|) , takes the form 

£d + 2??(V X ■ V y )G + (Hi ■ V x )v + V y x (v x G ) 

-(vV y )Ho = Ai(Ho + G ). (A.8) 

Substituting (|A.6| ) in ( |A.8| ) and averaging, we obtain 

V y x ( (v x S) H ) = AiH . (A.9) 
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So far, developments strictly followed those of Vishik (1987), but here 
the point of divergence is reached. Since v is parity-invariant (|5|) and S is 
parity-antiinvariant ( |A.7j ), we have (v x S) = 0, and therefore ( |A.9[ ) implies 
Ai = 0. Upon inserting ( |A.6| ) into the fluctuating part of ( |A.8| ), we obtain 

3 3/ aq \ pfVfk 

CG 1 = -(H x • V x )v + V V -2r,?P- + v m (S k + e k ) - vS ( m * ^° 



9x m J dy m 

(A.10) 



from which one finds 



3 3 3 q rrk 

G 1 = £ S*(x)ff*(y) + EE r mfe (x)^(y). (A.ll) 



k=l k=l m=l 



Here, the nine vector fields T mk are solutions to the second auxiliary equation 
(1131). Note r mfc (x) = — T mk (— x) by the symmetry properties of v, and of 
the operator C. By taking divergence of the second auxiliary equation one 
finds V x ■ T m fc + S™ = 0. Hence, for n = 1 (|A.3|) is automatically satisfied. 
in. The e 2 term of ( |A.4|) takes the form 



CG 2 + V (2(V X • V y )Gi + V 2 (H + Go)) + (H 2 • V x )v 
+ V y x (v x Gi) - (v • V y )Hi = A 2 (H + G ). (A.12) 
Upon substitution of (|A.11|) , averaging produces the second-order equation 

3 3 FlT-f^ 

AfHo = 7/V y Ho + Vy x E ( v x r^)/ = A 2 H , (A.13) 

k=l m=l ^ m 

which defines the eddy diffusivity operator. So, the leading terms of the ex- 
pansions of the mean magnetic field, Ho, and of the associated eigenvalue, 
A 2 , satisfy the eigenvalue equation for this operator. Its eigenvectors are 
Fourier harmonics: Ho = he iq y , where h and q are constant vectors satisfy- 
ing (P]) and flT5|). Note that the fluctuating part of the leading term of the 
eigenmode, G , is now completely determined by ( |A.6| ). 

At this point our construction goes beyond that of Lanotte et al. (1999) 
and, in the spirit of Vishik (1987), we show how to extend the analysis to 
higher orders in e. 

The fluctuating part of ( A.12 ) reads 

V \S J?i (^m dy m + ^ dxi dy m dyi) + V y G °) 
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:H 2 • V x )v + (v • V y )Hi 



/ 3 3 3 f)H k \ 

- V y x j> x S k )H* + E ( v x r m*> TT^ + A 2G . (A.14) 

\fc=l k=l m=l °^/ m / 

Here, only the terms involving H 2 and dH\jdy m are not yet determined. 
However, as one can immediately see, x-dependent prefactors in front of 
these terms are identical to those, which are in ( |A.10|) in front of Hi and 
dH^/di/m, respectively. This implies a representation 



3 3 3a TTk 

G 2 = ]T S fc (x)ff*(y) + EE r mfc (x)^(y) + Q 2 (x,y), 

k=l k=l m=l °Vm 

where the vector field Q 2 can be uniquely determined from 



/ 3 3 3 a-p a2 jrk N 

£Q 2 = - V 2 E E E + v y G 



mfc 

,,_ JW =i^ dx t dy m d yi 



3 3 FlT-f^ 

k =i m =i vy™ 

iv. As already stated, A 2 is one of the two eigenvalues of (p^)- ([T5|) . In 
what follows we assume that the other eigenvalue, A' 2 , is distinct from A 2 and 
we denote by h' the associated eigenvector. 

Equations arising at higher orders in e are now solved recursively. Assume 
that all the equations up to order e^ -1 have been solved, so that the following 
information has been obtained: 

• vector fields H n and G n for all n < N — 2; 

• representations of G n of the form 

3 3 3 a zrk 

G n = £ S fc (x)tf n fc (y) + E E r mfe (x)— ^i(y) + Q„(x,y) (A.15) 

k=l k=l m=l 

for n = iV — 1 and n = N — 2 with known vector fields Q n , (Q n ) = 0; 

• quantities A n for n < N. 

Upon substitution of ( |A.15|) for n = N — 1 the average of the equation at 
order becomes 

N-3 

(M - A 2 )Hat_ 2 - AtvH = -V y x (v x Qjv_i) + E A w _ m H m , 

m=l 
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where the right-hand side at this stage is a known vector field. Projecting 
this equation out in the direction of Ho one can uniquely determine A at. In 
the complementary invariant subspace the operator M. — A2 is invertible, 
and hence Hat_ 2 can be determined up to an arbitrary multiple of Ho- This 
additive term can be neglected: one can demand that Hat_ 2 belong to the 
complementary subspace. This is a normalization condition related to the 
fact that eigenvectors can be multiplied by arbitrary (analytic) functions of 
e. As a result, Gn-2 is now determined by (|A.15| ) with n = N — 2. 

Upon use of (|A.15|) for n = N — 1 the fluctuating part of the equation at 
order becomes 

£Q N — (2 V V ( dSfc d#Ar-i _j_ y> dr mk d 2 H%_ 2 \ 

dxi dyidy 



+2(V X ■ V y )Q A r_ 1 + V^G iY - 2 j - (Hjv • V x )v + (v • V y )Hjv_i 

/ 3 3 3 arrk \ 

"Vy X £( V X Sfc)^-! + £ E ( V X T mk} ~f^~~~ + i V X Q^-l} 
Vfc=l fc=lm=l / 

iV-2 

+ E Ajv_ m G m . (A. 16) 

m=0 

Only the terms involving H^v and derivatives of H^_ 1 are not yet determined. 
Like in ( |A.14| ), x-dependent prefactors in front of these terms are identical 



to those in front of Hi and derivatives of Hq in (|A.10|) , respectively. Hence, 
( |A.16[) yields a representation of the form (|A.15| ) for Gat. An equation for 



Qn can be obtained from ( |A.16|) by replacing by Qat and omitting all 
terms involving H^ and derivatives of H^-i- 

Thus, a complete formal decomposition of magnetic modes and their 
eigenvalues has been constructed. A step-by-step analysis of the solution 
reveals that 

H = he iq ' y , H n = XnhV qy Vn>0 
(with suitable scalars Xn) an d 

G n = g„(x)e 4q - y Vn>0; (A.17) 

thus the expanded eigenmode admits a representation (JJ) in the form of a 
plane wave in the slow variable. (For this reason it was sufficient to demand, 
for constructions presented under heading iv, that the two eigenvalues A2 
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and A' 2 of M. be distinct in the subspace of solenoidal vector fields ce* qy with 
c =const, instead of requiring that M. — A2 be invertible in the whole domain.) 
The plane-wave representation does not come as a surprise, since the domain 
of the magnetic induction operator M. splits into invariant subspaces, each 
comprised of vector fields of the form (|J) and categorized by wavevectors q. 

If the homogenized operator M. turns out to be elliptic, one can prove, 
following Vishik (1986, 1987), that any point A 2 of the spectrum of M. is 
associated with an analytic branch (|8]) of eigenvalues of the original magnetic 
induction operator A4. The proof relies on the general perturbation theory 
for linear operators (Kato, 1980); it will not be presented here. 

B Strongly negative eddy diffusivities 

In this Appendix we discuss a mechanism for appearance of negative eddy 
diffusivities with arbitrarily large magnitudes. 

Let i p (rj) denote the basis of small-scale 27r-periodic magnetic modes with 
vanishing mean and Cp(v) ~ the associated eigenvalues: 



(assuming in this section, for the sake of simplicity, that C is diagonalizable). 
Suppose that when molecular diffusivity rj is decreased, small-scale magnetic 
fields with vanishing mean lose stability in such a way that the dominant 
eigenvalue of the magnetic induction operator C passes through zero. Let r\ c 
be the critical diffusivity: 



Generically, the loss of stability occurs either in the proper subspace com- 
prised of parity-invariant fields, or in that of parity-antiinvariant fields, im- 
plying two possible variants of the mechanism. 

Suppose small-scale parity-invariant magnetic fields are the first to lose 
stability. Solutions S (77) to the first auxiliary problem ([12]) continuously 
depend on rj for 77 < rj c . The second auxiliary problem (|T3|) can be reexpressed 



■v 



fa)| = l, V-f p (r/) = 



Re C P (v) < Vp, V?? < Vc - Cifoc) = 0. 



as 



(B.l) 



p p 
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where 



"Ml 



and 

r mfc (77) = ^7mfc,p(?7)fp(?7) 

p 

are decompositions of the r.h.s. of ( p~3| ) and of r mfc (?7) in the basis of magnetic 
modes; summation over all parity-invariant eigenvectors is assumed in the 
two sums. From ( |B.1| ), •y m k, P (v) = K mk, P (v) / Cpiv) ■ Since generically the 
coefficients n m k,i{Vc) do not vanish, we have in the limit rj —*■ rj c 

Then the equation (|14D for the eigenvalue A2 asymptotically reduces to 

^|EE 'Wfac^V" = -A 2 h, (B.2) 

fe=lm=l 

where A(q) = q x (v x fi (?7c) ) ■ Equations ( |B.2|) and ( |I5|) yield two eigen- 
values. It is easily checked that one of them, denoted A' 2 , vanishes and thus 
is of no interest for us; the associated eigenvector h' is orthogonal to the 
wavevector q and satisfies 

3 3 

E £ ^mk,i(Vc)h' k q m = 0. 

k=l m=l 

The second eigenvalue is 

A 2 = o(q) = EE ^,i(^)^(q)g m , 

the associated eigenvector being h = A(q) (we consider a generic case, where 
A does not vanish). Therefore, if a = max| q | =1 a(q) > 0, then, in the limit 
77 — ► r/ c , the minimum magnetic eddy diffusivity becomes arbitrarily large 
negative: 



32 



Alternatively, suppose now that the loss of stability occurs in the subspace 
of small-scale parity-antiinvariant magnetic fields. Then the mechanism is 
similar. In this case already the solutions S(r^) to the first auxiliary problem 
( |12D tend to infinity. Indeed, flip ) can be represented as 



where 



and 



■^ = EM f i) f P (')) 



V 



are decompositions of the r.h.s. of the first auxiliary problem (|T|) and of 
Sk{i]) in the basis of magnetic modes; summation over parity-antiinvariant 
eigenvectors is assumed in both sums. Thus, Sk, p (rj) = Pk, P {v) / (p{v) ■ Since 
generically Pk,i(Vc) 7^ 0, we have in the limit rj — ► rj c 

Let 7 m denote the solution to 

ox m 

Then ^ 



and (|Tj) asymptotically reduces to 



B(q) EPM(^c)/»* = -A a h, (B.3) 



where 

B(q)= Eg m qx (vx 7m ). 



m=l 
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Equations ( |B.3| ) and (p~5| ) yield two eigenvalues. Again, one can easily check 
that one of them, A' 2 , vanishes; the associated eigenvector h' is orthogonal to 
q and satisfies 



k=l 



The second is 



fc=i 



with the associated vector h = B(q) (in the generic case B^O). Therefore, 
if b = max| q i = i fe(q) > 0, when 77 — > r] c the minimum magnetic eddy diffusivity 
becomes again arbitrarily large negative: 



Redely ~ -T-TT ^ -OO ■ 

C Chebyshev iterative methods for numeri- 
cal solution of systems of equations 

We discuss here three Chebyshev iterative methods suitable for large sys- 
tems of linear or non-linear equations. They employ extremal properties of 
Chebyshev polynomials and their roots. The first method (Chi), discussed in 
Subsection |C.1|, is applicable to problems with negative-definite self-adjoint 



Jacobians. An extension (Chi a), discussed in Subsection |C.2| , is an upgrade 



of Chi, which eliminates these restrictions. The second method (Ch2), dis- 
cussed in Subsection |U.3| , is designed for discretizations of elliptic differential 
equations, such as the auxiliary problems ([12]) and fll3|) . Chi and Chi a were 
introduced by Podvigina and Zheligovsky (1997) and Ch2 is presented in 
detail in Zheligovsky (2001). For general background on Chebyshev iterative 
methods, see Axelsson (1996). 

C.l The basic method: restricted case {Chi) 

We consider numerical solution of a system of equations 

F(z) = 0, (C.l) 
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where F : — > C N is smooth and the number of equations, N, is large. 

Let z = Z be a solution to ( |C.1| ). We first consider the case, where 
the Jacobian matrix F'(Z) is diagonalizable and its spectrum is (real and) 
negative. Then the solution can be computed by iterations 

zi = z + /iF(z ), (C.2) 

Zfc+i = lk (z fc + h¥(z k )) + (1 - 7fc)z fc _i, k > 0, (C.3) 
provided z is chosen in the basin of attraction of Z. Here 

9 _ 9 TM _ A + A 

7fc = 2/i- — — , \x = -, (C.4) 



A + A' ,K ~ ^T h+1 (fiy ^"A-A ; 

where < A < A for any A > |F'(Z)| (the optimal value A = |F'(Z)| is 
usually unknown), and where the TVs are Chebyshev polynomials of the 
first kind. The 7fc's can be evaluated recursively: 

Consider how the discrepancy vectors £ fc = — Z evolve in the course of 
iterations ( |C.2| )-( |CT5|) . Introduce the eigenvalues rji of F'(Z) (here, 
> rji > —A) and the associated eigenvectors {ej}, and decompose the 
initial discrepancy in this basis: 

N 

Co = E c ' e - 

i=l 

It is then easily shown that 

^ T fc (/i(l + %)) ■ 
i=i Jfcl^J 

(the o(]Co I) correction does not appear if the mapping F is linear). 

At this point, the standard strategy is to use the least deviation principle 
for Chebyshev polynomials: the iterative process (|C.2|) -( |C75| ) with A > 



ensures maximum uniform decay of the discrepancy vector component in 
the proper subspace of F'(Z) associated with the part of the spectrum in 
the interval —X>f]i> —A. The alternative procedure used by Podvigina 
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and Zheligovsky (1997) relies, in addition, on the extremal property of the 
largest root of a Chebyshev polynomial: for A = 0, the iterative process 
gives optimal reduction of the discrepancy vector component in the proper 
subspace of F'(Z) associated with eigenvalues rji closest to zero. Therefore, 
a suitable algorithm, called here Chi, consists of a succession of iterative 
sequences flC3D- (|d5D with A > and A = 0. 

As usually, lowering of the condition number of the system allows to 
improve efficiency of the algorithm. If, in addition, F'(Z) is self-adjoint, such 
lowering can be achieved by replacing (|C.1| ) with an equivalent system 

F(z) = 0, 

where 

F(z)=P(z).F(z), (C.6) 

P(z) being a self-adjoint and positive-definite linear preconditioning opera- 
tor. It may be checked that for such a P(z) the mapping F has the required 
properties stated at the beginning of this subsection. 

C.2 An extension of Chi: general case (Chi a) 

In the case of an arbitrary smooth mapping F, to which we turn now, a 
modified system 

F(z) = 

can be considered, where 

F(z) = -P 1 (z).(F'(z))*-P 2 (z)-F(z), (C.7) 

Pj (i = 1,2) being arbitrary positive definite self-adjoint operators and 
M* = M* being the operator adjoint to M. The method discussed above 
can be applied for the solution of the modified system, since F'(Z) has all 
the required properties. 

Termination of individual sequences and selection of values of A in suc- 
cessive sequences is controlled by the following rules: 

1. In the first iterative sequence A = Ai. 

2. In any of the following iterative sequences A = 0, unless A was zero in the 
orevious sequence and that sequence was short (i.e. it consisted of less than 

71 \J sXT ^ derations, where [•] denotes the integer part of a number), in 
which case in the next sequence A = Ai. 
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3. An iterative sequence with A = Ai can be terminated only after a pre- 
scribed number of iterations have been performed. 

4. In an iterative sequence with A = the errors |F(z^)| can initially grow 
monotonically, if the respective errors |F(z^)| decrease monotonically. (An 
increase of |F(zfc)| at the first iteration indicates that A is underestimated.) 

5. After the initial phase allowed by rules 3 or 4 is completed, an iterative 
sequence is terminated if both discrepancy reduction rate values, | In p^°] 

and | In ^ i are smaller than at the previous iteration. 

6. Either the last or the next to last iterate is accepted as the output of a 
sequence - specifically, the one, for which the error |F(zjk)| is smaller. 

With these rules, the behavior of the algorithm proved insensitive to the 
particular choice of Ai in the interval 0.1A < Ai < 0.5A. 

Only trivial changes to this set of rules are required to make it applicable 
to Chi. 

Chla proved robust and efficient, e.g. in computation of steady states of 
the Navier-Stokes equation (Podvigina, 1999a,b). However, the use of ( |C.7|) 
may be potentially disadvantageous, because: 

• The presence of the adjoint to F'(z) leads to an increase in the condition 
number, which may not be completely neutralized by the preconditioning 
operators Pi and P 2 . 

• Computation of F involves more operations (and hence larger codes) than 
that of the original mapping F. 

• Spurious roots may appear. 



C.3 An adaptation of Chi to elliptic PDE's (Ch2) 

We present here an alternative to Chla, better suited to problems arising 
from discretization of elliptic partial differential equations, when they involve 
not just leading-order elliptic operators (such as the Laplacian), but also 
lower-order non-self-adjoint perturbations. We do not attempt to present 
this method in full generality (see Zheligovsky, 2001). In what follows the 
resolution of the discretization, assumed to be large, is the control parameter, 
which we denote by R. 

Construction of the modified mapping ( |U.7| ) can be avoided under the 
conditions, that the basic mapping F can be split: 

F(z) =D(z) + A(z) (C.8) 
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and if a linear invertible preconditioning operator P can be selected, such 
that: 

1. All eigenvalues cXj of the linear operator T> = P • D'(Z) are negative. 
(If P is a self-adjoint positive-definite operator and D'(Z) is a self-adjoint 
negative-definite operator, this condition is automatically verified.) 

2. The ratio 5 = \A\/A is small for large R. Here A = P ■ A'(Z), and 
A is a "tight" upper bound for the spectral radius of V, i.e. a bound, which 
stays of the order of the spectral radius for large resolutions. 

3. The operators T> and A nearly commute: \T>A — AT>\ = o(\VA\). 
(Note, that this condition is not that restrictive: it holds automatically in 
the absence of preconditioning, if T> is for example the Laplacian.) 

A decomposition ( |C.8| ) arises naturally if the problem ( |C.1| ) is a discretiza- 
tion of a dissipative partial differential equation: then the mapping D can 
be identified with the higher-order elliptic part of the equation, and A - 
with the remaining terms. Auxiliary problems flT2| ) and (|H|) arising from the 
two-scale expansion fall into this class. 

Under the conditions stated above a modification of the method discussed 
in Subsection |C.1| , named Ch2, can be applied. Solution proceeds by a suc- 
cession of Chebyshev iterative sequences ( |CJ.3| )- (|^T75D , but the first iteration 
in all sequences is made differently. Suppose a termination condition (for 
example, the one considered in the previous subsection) for a Chebyshev se- 
quence number n — 1 is verified for the first time in this sequence for the 
iterate z^- -1 (we then call K the termination signal number). Then the se- 
quence is continued by M further iterations, and we start the next sequence, 
by selecting 

z£ = z- 1 ; z? = z^ 1 + h (V(zl~ l ) + A(*%2«)) ■ 
Subsequent iterations are defined as in flC.3| ): 

z£ +1 = 7fc (z™ + hF(z n k )) + (1 - 7fe )z]U k > 0. 

(The quantities h, \i and 7^ were defined by ( |C.4|) - (|C.5|) and F - by (|C.6| ).) 

In this iterative process, errors behave in the following way. Discrepancy 
vectors Ck = z k ~ ^ can be represented in terms of Chebyshev polynomials 
of the first (Tfc) and of the second ([/&) kind: 

C = ^(UQ)a n + U k (Q)b n ) + o(|C|), 
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where Q = + h(V + A)), I is the identity operator, and vectors a n and 
b n , which we call the basic error vectors, can be determined from £q and Ci ■ 
(The terms o(|Cfc|) do not emerge if the mapping F is linear; in what follows 
they are neglected.) 

Now, suppose that for the (n — 2)-nd sequence the termination signal 
number K was the same as for the (n — l)-st iteration (in computations, 
successive Chebyshev sequences often do have equal lengths). Then it can 
be shown that 

• |b n | = 0(5|a n_1 |) and thus convergence of b n is slaved to convergence 
of a n . 

• Within the succession of iterative sequences, the basic error vectors 
evolve as follows: 

N n-1 

a " = E ^VtCM^ + h^T^Ae,) + o^a"" 1 !), (C.9) 
i=i 

where 

f = 1 + 1 ( T ^+Mfe) _ .\ 

£i = /i(l + h(Ti), Oi being eigenvalues of V, and the coefficients a™ -1 can be 
found from the decomposition of a n_1 in the basis of eigenvectors ej of T>\ 

N 

a**" 1 = Ear 1 *- 

i=i 

Thus, reduction of the basic error vectors a n is predominantly controlled 
by V (note that h\A\ = 0(5)). The standard first Chebyshev iteration ( |C.2| ) 
is recovered when one sets M = and hence = 1. It may be shown that 
our procedure leads to factors < fi < 1, and thus to enhanced convergence, 
for eigenvalues sufficiently small in absolute value (provided A = or A is 
small). For the other eigenvalues good convergence is ensured as usually. 
(Actually, because of the presence of the non-self-adjoint perturbation A, 
fast convergence is guaranteed only for eigenvalues large in absolute value; 
for those of intermediate magnitude slowdown might occur, but has not been 
observed in any of the cases we studied.) Of course, these convergence ar- 
guments hold only in so far as the resolution R is high enough and thus 5 is 
small enough to justify keeping only the leading terms in ( |C.9] ). 
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In the application of the generalized Chebyshev method to the auxil- 
iary problems dl2|) and ([13|) , we take D = 77 V 2 (more precisely, D is the 
Fourier-Galerkin discretization of this operator). We select the precondi- 
tioning operator P = (— D)~ 9 for some < q < 1 by analogy with Podvigina 
and Zheligovsky (1997) and Podvigina (1999a,b). Then V = -(-D) 1 "" 
is a pseudodifferential operator of order 2 — 2q, and A - of order 1 — 2q. 
Hence, Condition 2 is inherited from the non-discretized versions of the op- 
erators for q < 1. Condition 3 is satisfied for q < 3/4. Indeed, in the 
commutator T>A — AT> the leading symbols cancel out and therefore its or- 
der is by 1 less than the order of each of the operators T>A and AD; thus 
\VA-AV\ = 0( J R max (°' 2 - 4( ')), \VA\ = O(R m!lK(0 ' 3 -^), where R is the maxi- 
mum wavenumber of the Fourier harmonics, retained in the abridged series, 
this completing the argument. Since Condition 3 is also trivially satisfied for 
q = 1, it is natural to consider the larger interval < q < 1, and we find the 
convergence to be the best for q = 3/4. When applied to (|i~2| ) and (|T3|), Ch2 
proved superior to Chi a. 
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